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We report small-angle neutron scattering measurements of the vortex lattice (VL) structure in single crystals of the 
lightly underdoped cuprate superconductor YBa2Cu306.85- At 2 K, and for fields of up to 16 T applied parallel to the 
crystal c-axis, we observe a sequence of field-driven and first-order transitions between different VL structures. By 
rotating the field away from the c-axis, we observe each structure transition to shift to either higher or lower field 
dependent on whether the field is rotated towards the [100] or [010] direction. We use this latter observation to argue 
that the Fermi surface morphology must play a key role in the mechanisms that drive the VL structure transitions. 
Furthermore, we show this interpretation is compatible with analogous results obtained previously on lightly overdoped 
YBa 2 Cu 307 . In that material, it has long-been suggested that the high field VL structure transition is driven by the nodal 
gap anisotropy. In contrast, the results and discussion presented here bring into question the role, if any, of a nodal gap 
anisotropy on the VL structure transitions in both YBa2Cu306.85 and YBa 2 Cu 307 . 


1. Introduction 

It is now more than fifty years since the seminal theoret¬ 
ical work of Abrikosov wherein the superconducting vortex 
lattice (VL) was first shown to be expected in type-II super¬ 
conductors.* This pioneering work sparked the intense and 
continuing research efforts into the physics of vortex matter, 
which assumed even greater importance with the discoveries 
of High-Tc cuprate and other unconventional type-II materi¬ 
als with novel pairing symmetries. A full understanding of 
the properties of vortex matter is also crucial from the techno¬ 
logical perspective because type-II superconducting materi¬ 
als are, for example, key constituents of high current-carrying 
applications. From the latter viewpoint, it becomes vital to 
understand the vortex properties of High-Tc cuprates that dis¬ 
play superconducting critical temperatures above those of the 
77 K boiling point of liquid nitrogen, since this is a much 
more economical cryogen than liquid helium. 

Close to optimally-doped YBa2Cu307_5 is a prototypical 
High-Tc cuprate with superconducting critical temperatures in 
the region of 90 K, the exact value being composition depen¬ 
dent.^ Indeed, the VL properties of YBa2Cu307_5 have been 
the subject of numerous investigations over the years,with 
the bulk probe of small-angle neutron scattering (SANS) play¬ 
ing an important role in revealing the effects of effective mass 
anisotropy, VL melting and pinning, and field-induced non¬ 


locality.' Despite the large body of work, it is only in 
the last decade that SANS observations of the intrinsic mag¬ 
netic field- and temperature-dependent VL structures could be 
achieved.The more recent achievements can be largely 
attributed to significant improvements in sample quality, since 
these observations were not complicated by twin-plane pin¬ 
ning effects which hampered the earlier efforts to unveil the 
intrinsic VL properties.*' *^’ *^ 

In particular, at low temperature in detwinned and lightly 
overdoped YBa2Cu307, the field-evolution of the VL struc¬ 
ture was recently studied for magnetic fields up to 16.7 T ap¬ 
plied parallel to the crystal c-axis.^*“^^ The most striking ob¬ 
servations include successive field-driven first-order reorien¬ 
tation transitions of the VL coordination at both 2.3(2) T and 
6.7(2) T,^* and a smooth evolution of the high-held VL that 
passes through the perfectly square VL coordination expected 
to be stabilised by a pure (7-wave gap symmetry. In general, 
anisotropies in both the Fermi surface and the superconduct¬ 
ing gap can be responsible for the field-dependence of the 
precise VL structure,and also cause VL structure tran¬ 
sitions. However, and as discussed in Refs.,^^’^** YBa2Cu307 
remains a high-profile example of an unconventional ^7-wave 
superconductor in which the physical origins of the VL struc¬ 
ture transitions are not completely understood. 

Towards addressing this problem, here we present 
new SANS studies of detwinned and lightly underdoped 
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YBa2Cu306.s5- Despite the light underdoping which intro¬ 
duces oxygen vacancies, and hence pinning centers, into the 
CuO chains, we nonetheless observe well-ordered VLs up to 
the highest magnetic field of 16 T applied parallel to the c- 
axis. We observe a sequence of VL structure transitions in this 
material that is analogous to that seen in YBa 2 Cu 307 , though 
the transition fields are slightly modified. We also present ex¬ 
tensive measurements of the VLs observed in YBa2Cu306.85 
as a function of angle of applied field to the c-axis. The VL 
structure transition fields are found to display a strong and 
characteristic angle-dependence, and we use this observation 
to provide new insights into the physical origin of the two 
field-driven VL structure transitions in YBa 2 Cu 307 -( 5 . 

2. Experimental Method 

Single crystals of YBa 2 Cu 306.85 were prepared by the 
flux growth method, and from a starting mixture of BaC 03 , 
CuO and Y 2 O 3 powders molten in yttrium-stabilized Zr 02 
crucibles.The as-grown crystals are crystallographic ally 
twinned and also have an ill-defined oxygen content. The in¬ 
dividual crystals were therefore detwinned at elevated tem¬ 
perature using a uniaxial stress applied along a (100) axis. To 
obtain the YBa 2 Cu 306.85 composition, the crystals were sub¬ 
sequently annealed at 520° C in an O 2 atmosphere. 

The sample used for the SANS measurements was a co¬ 
aligned mosaic of eight YBa 2 Cu 306.85 single crystals of total 
mass 85 mg. Using single crystal neutron diffraction to study 
the (200) and (020) nuclear reflections, the bulk twin-domain 
population ratio was determined for each individual crystal. 
The YBa2Cu306.85 single crystals were found to be almost 
completely twin-free, and the entire mosaic displayed a de- 
twinning ratio of 98.6(5)%. The superconducting critical tem¬ 
perature (Tc) of each individual crystal was also measured us¬ 
ing a Quantum Design physical property measurement system 
(PPMS). The overall sample displayed a Tc of 90.5(1.0) K, 
where the error represents the spread in amongst the crys¬ 
tals. Comparison between this Tc value and the data reported 
in Ref.^ gives a sample doping level of b = 0.15(1). The sam¬ 
ple mosaic was co-aligned and mounted on a 1 mm thick Si 
plate with the crystal c-axis of the sample perpendicular to the 
plate. 

The SANS experiments were conducted on the SANS-I 
instrument at the Swiss spallation neutron source (SINQ), 
Paul Scherrer Institut (PSI), Villigen, Switzerland, the D11 at 
the Institut Laue-Langevin (ILL), Grenoble, France, and the 
NG3-SANS instrument at NIST Center for Neutron Research, 
USA. The neutrons were collimated over 6-18 m (SANS-I), 
5.5 m (Dll) and 4.5-8m (NG3-SANS), and wavelengths be¬ 
tween 4.5-10 A were selected with a typical FWHM spread 
of 10%. A position-sensitive 2D multi-detector placed be¬ 
hind the sample was used to detect the scattered neutrons. 
At SANS-I, Dll, and NG3-SANS, 11 T, 17 T and 9 T hor¬ 
izontal field SANS cryomagnets were respectively available. 
In each cryomagnet the sample was initially installed with the 
crystal c-axis parallel to the applied field direction, and ap¬ 
proximately parallel to the neutron beam. For the 11 T mag¬ 
net at PSI, an additional motor-head allowed rotation of the 
sample stick about the vertical axis inside the variable tem¬ 
perature insert (VTI) of the cryomagnet. This motor provided 
in-situ control of the direction of the magnetic field with re¬ 


angle-dependent measurements of the VL structure presented 
in Section 3.2. 

At each SANS instrument the cryomagnet was placed on 
a goniometer table that enabled both rotation and tilt of the 
sample-magnet ensemble. The SANS measurements of the 
VL were done by rotating the sample and cryomagnet to¬ 
gether through a range of angles that moved the VL diffrac¬ 
tion peaks through the Bragg condition. By summing over all 
measurements at different angles, the SANS diffraction pat¬ 
terns shown in this paper are obtained, and wherein all the 
VL diffraction spots can be observed in a single image. 

For all the low temperature measurements reported in this 
paper, the VL was prepared using the oscillation-field-cooling 
(OFC) method. By this method the sample is field-cooled 
through Tc until the measurement temperature in a time- 
dependent magnetic field that weakly oscillates around the 
target field. The oscillation amplitude used was typically 1 % 
of the target field. The OFC method was used to prepare the 
VL since previous work on twin-free YBa 2 Cu 307 showed 
such an approach allows i) the VL to achieve an improved 
perfection, and ii) a preferred VL coordination due to a lower 
effective temperature than achieved when cooling in a static 
field.Additional SANS measurements not reported here 
confirmed that the OFC method produced a similar effect on 
the VLs observed in the present YBa2Cu306.85 sample. Back¬ 
ground SANS measurements were conducted at T = 95 K 
and subtracted from the low temperature measurements in or¬ 
der to leave just the VL signal. SANS data visualization and 
analysis was carried out using the GRASP software.^® 

3. Results 

3.1 Vortex lattice structure for hqH || c in YBa2Cuj,0(,,is 

Figs. l(a)-(c) show a selection of VL SANS diffraction pat¬ 
terns obtained from YBa2Cu306.85 at T = 2 K, and at different 
magnetic fields applied parallel to the crystal c-axis {} 1 qH || c). 
Fig. 1(a) shows the diffraction pattern obtained at 1 T that is 
typical for fields up to 2.37(2) T. As shown schematically in 
Fig. 1(e), this VL is composed of two distorted hexagonal do¬ 
mains that each have two spots aligned with a {110} plane. 
The alignment of VL domains along these directions is a clear 
indication that the VL orientation is not intrinsic, but instead 
controlled by the small volume of twin-planes remaining in 
the sample. Consequently, we refer to this two domain 
VL structure as structure TP. As illustrated for one of the two 
domains in Fig. 1(e), the primitive cell opening angle (f>\ is 
almost bisected by the horizontal b*-axis, but not exactly. 

In the field range between 2.37(2) T and 5.37(2) T, structure 
TP coexists with an additional single distorted hexagonal VL 
domain, which we call structure B. Fig. 1(b) shows the typical 
SANS pattern for this field range, with a schematic illustration 
of all the VL domains present shown in Fig. 1(f). In contrast 
to structure TP, structure B is aligned intrinsically with the 
crystal axes, and its primitive unit cell is bisected exactly by 
the vertical a*-axis. Across this field range the VL domains of 
both structures TP and B are always observed to co-exist and 
neither exclusively populates the entire sample. By surveying 
the individual crystals of the sample mosaic, we confirmed 
the structure co-existence in all crystals, though the relative 
populations of each structure differed slightly. This observed 
variation could reflect small differences in either oxygen stoi- 
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Fig. 1. (Color online) SANS diffraction patterns obtained from the VL at 
2 K and (a) || c = 1 T, (b) iiqH || c = 4 T, (c) fi^H || c =16 T, 

and (d) =1 T applied 10° away from c. Panels (e) to (h) show schematic 

drawings of the results shown in (a) to (d), and which clarify the different VL 
domains present. Panels (e) and (f) each show VL domains, denoted TPl and 
TP2, whose alignment is affected by pinning effects due to the small residual 
density of twin-planes. In contrast, structures B (panel (f)) and A (panels (g) 
and (h)) are aligned with the crystal axes and their alignment is considered to 
be intrinsic. 


spacings between the residual twin planes. 

Increasing the field beyond 5.37(2) T leads to a suppres¬ 
sion of both structures B and TP. Instead the VL structure 
becomes composed of a different intrinsic and single VL do¬ 
main aligned with the crystal axes that we call structure A. 
Structure A prevails until the highest measured field of 16 T, 
the SANS diffraction pattern of which is shown in Fig. 1(c). 
This structure is sketched in Fig. 1(g), and we define the prim¬ 
itive unit cell to be bisected exactly by the horizontal b*-axis. 
A general feature of the data obtained in the high-held range 
is that the high-held diffraction pattern is dominated by a 
rhombic VL coordination characterized by four intense Bragg 
spots, in contrast to the hexagonal VL domains observed at 
lower helds. 

Fig. 2 shows the //o77-dependence at T = 2 K of the prim¬ 
itive cell opening angles for structures TP, B, and A. For the 
TP structure, (ft is an average taken over the two constituent 
domains. Overall, with increasing /tq// there is a smooth evo- 



Fig. 2. (Color online) The VL opening angles of the various VL primitive 
cells for fifiH || c and at T = 2 K. The dashed line is a guide for the eye. 
The shaded region denotes the field-range of co-existence for structures TP 
and B. Inset: the iii^H || c-T VL structure phase diagram. The symbols denote 
the transition fields between the different VL structure phases, as determined 
by tracking the VL structures in the vicinities of the phase boundaries. The 
melting line was obtained from Ref..^** 


order re-orientation transition which occurs at 5.37(2) T into 
the high-held phase where only structure A exists. At higher 
held, the primitive cell of structure A evolves smoothly with 
increasing up to 16 T, and at 10.6(2) T it passes through 
the value of 90° expected for a perfectly square coordination. 
The co-existence of structures TP and B between 2.37(2) T to 
5.37(2) T is also seen clearly, while for helds outside this held 
range just a single VL structure type populates the sample. 

Bearing in mind the similarity between the primitive cell 
orientations for both structures TP and A, structure A can be 
expected to be the intrinsic structure that would be present 
at lower helds in the absence of an orientational effect due 
to twin-plane pinning. This assumption can be readily tested, 
since it is well-known that in lightly-twinned YBa 2 Cu 307-5 
samples, twin-plane pinning effects can be suppressed by ro¬ 
tating the sample so that all twin-planes form a small angle of 
~5-10° with respect to the held direction.After rotating 
the present YBa 2 Cu 306.85 sample by 10° around the vertical 
a*-axis with respect to a 1 T held, we hnd that such a rotation 
is indeed sufficient for the VL to overcome the twin-plane pin¬ 
ning. This is seen in Fig. 1(d) which shows that in the rotated 
condition the VL diffraction pattern is composed of just a sin¬ 
gle distorted hexagonal domain aligned with the crystal axes, 
and not the {110} planes. As illustrated in Fig. 1(h), the prim¬ 
itive cell for this VL structure can be treated in the same way 
as structure A, with the denoted primitive unit cell bisected by 
the horizontal axis. 

Returning to fioH || c, another important feature seen in 
Fig. 2 is the broad held co-existence range of the TP and B 
structures. Such a co-existence exemplihes the near degen¬ 
eracy between the two structure-types, and also shows that 
at 2.37(2) T a hrst-order-type of VL structure transition oc¬ 
curs between a fraction of the TP structure and structure B. 
Since such a ‘partial’ transition occurs within all crystals of 
the sample mosaic, it could be argued that the particular oxy¬ 
gen doping displayed by the sample appears to be close to a 
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Fig. 3. (Color online) Sketches to illustrate the expected effect on the ob¬ 
served VL structures for SANS measurements done with the field applied 
at an angle to the crystal c-axis. In (a) the sample is rotated relative to the 
^oH axis by an angle a ai'ound a*. In (b) the sample is rotated by an angle 
around b*. On the right side of each panel are sketches of the expected SANS 
VL diffraction patterns for the different cr/yS, and where directions qx and qy 
respectively denote the horizontal and vertical directions on the detector. The 
filled circles denote the Bragg peaks due to structure A, while empty circles 
those of structure B. 


critical level for the occurrence of a complete phase transition. 
We also point out that this coexistence phase appears not to be 
a property of twin-plane pinning. As will be shown later, by 
again rotating the sample with respect to the field direction to 
suppress the pinning, the broad co-existence region survives, 
and is between structures A and B instead of structures TP 
and B. 

The inset of Fig. 2 shows the ^qH || c-temperature VL 
structure phase diagram deduced from measurements where 
we tracked the VL transition fields. Here it is seen that 
the VL structure phase boundaries are only very weakly T- 
dependent, and so the field-width of the co-existence region 
remains essentially unchanged until the melting temperature. 
This behavior is similar to that reported for YBa 2 Cu 307 in 
Ref.^^ with the same discussion presented there applicable 
also here for YBa 2 Cu 306 . 85 . 

3.2 Vortex lattice structure for /jqH not parallel to c in 

YBa2CuiO(,.s5 

To provide further insight into the //o^f-dependent be¬ 
haviour of the VL observed for poH || c, we conducted mea¬ 
surements of the VL structure for a range of angles between 
the applied field and the crystal c-axis. As illustrated in Fig. 3, 
we explored two directions of sample rotation relative to poH: 
a) rotating the field by an angle a around the a*-axis (towards 
[010]), and b) a rotation by an angle (3 around the b*-axis (to¬ 
wards [100]). 

In Fig. 3 we also illustrate the change in distortion of hexag¬ 
onal VL structures expected from anisotropic London theory 
as a consequence of the sample rotation. For pqH || c, i.e. 
a = y 6 = 0° we have seen in Sec. 3.1 that, although the orien¬ 
tation is determined by pinning, the low field hexagonal VLs 
are distorted along the a*-axis away from an isotropic hexag¬ 
onal coordination. Anisotropic London theory provides the 
simplest interpretation of the distortion in terms of an effec¬ 
tive mass anisotropy since, 77 = jab — 3alAb oc 
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Fig. 4. (Color online) The VL structures observed at T = 2 K, and for 
fiQH = 4T for (a) a = 30°, (b) o' = 10° and (c) yS = 20°. 


Bragg spots in a hexagonal VL domain. Since rj is related ge¬ 
ometrically to the VL opening angles defined in Fig. 1, it is 
readily seen that 77 > 1 at low field where London theory is 
most valid. In the zero-field limit we find 77 = 1.25(3) for 
YBa 2 Cu 306 85 ^ which indicates a lower effective mass, and 
hence higher super-current density, along the b*-axis com¬ 
pared with the a*-axis; an observation consistent with the ex¬ 
pected existence of superconducting CuO chain states. 

This suggestion is further supported by noting that the value 
of j] - 1.25(3) for YBa 2 Cu 306 . 85 , which has oxygen vacan¬ 
cies in the CuO chains, is less than that of 77 = 1.29(2) found 
for YBa 2 Cu 307 which has completely filled CuO chains,^'' 
and significantly less than the value of 77 = 2.57(5) reported 
for the stoichiometric double-chain compound YBa 2 Cu 408 .^^ 
By increasing either a or /3 the supercurrent response in the 
plane orthogonal to the field includes a contribution due to m*. 
Since 777 *» m* and m^, increasing either a or/3 always leads 
to an increase of the effective mass in the horizontal direc¬ 
tion. Therefore, and as shown in Fig. 3, increasing a causes 
the hexagonal VL distortion to decrease, and so a reduction 
of angle f characteristic of structure A, or an increase in an¬ 
gle p characteristic of structure B. The converse is realised for 
increasing f since the VL distortion will increase, with f in¬ 
creasing and p decreasing. It should be noted that even while 
anisotropic London theory provides a handy guide for under¬ 
standing the angle-dependence of the VL distortion, no field- 
driven VL structure transitions are expected within the theory 
for fields applied away from perpendicular to the basal plane. 
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Fig. 5. (Color online) Magnetic field-Of/yfi VL structure phase diagram at 
7 = 2 K. Each symbol denotes where the VL structure was measured by 
SANS. The phase boundaiy lines are estimated from the data. The dotted 
line close to zero angle {p-oH || c) and for fields below 5 T indicates the 
approximate region where vortex pinning can determine the VL orientation. 


Fig. 6. (Color online) The /io^f-dependence of the VL opening angle of (a) 
stmcture A and (b) stiaicture B over a range of a and /5 rotation angles. All 
data were obtained at T =2 K, and all dashed lines are guides for the eye. The 
legend shown in (b) applies to both panels. 


be understood as a consequence of higher-order anisotropies 
beyond those included in local anisotropic London theory. In 
what follows, it is observed that the //o^f-dependence of VL 
structure transitions are indeed strongly dependent on the di¬ 
rection of crystal rotation with respect to the field. 

The usefulness of VL measurements with hqH at an an¬ 
gle to the crystal c-axis is exemplified by the data shown in 
Fig. 4. Here we present VL diffraction patterns all obtained 
in fioH=4 T but for different values of a and (3. Beginning at 
II c. Fig. 1(b) showed the VL at 4 T to be composed of 
structures TP and B. In contrast, with a - 10°, an angle where 
pinning to twin-planes is suppressed, we observe the coexis¬ 
tence of structures A and B [Fig. 4(b)]. Further increasing a 
to 30° [Fig. 4(a)] leads to a suppression of structure A, and 
structure B becomes the only VL structure type in the sam¬ 
ple. On the p side the inverse tendency with angle is observed. 
Fig. 4(c) shows that for p - 20°, structure B is entirely sup¬ 
pressed, and structure A is the only structure-type in the sam¬ 
ple. Therefore, at a constant field of 4 T we observe an angle- 
driven first-order VL structure transition between structures A 
and B. From the most general perspective, this observation is 
in agreement with the expectations of the ‘hairy-balT theorem 
as applied to VLs, and which states that at any field discontin¬ 
uous VL transitions must exist as the angle of the field to the 
crystal is varied.^®’® 

In Fig. 5 we show an overall angle-/ro^f VL structure phase 
diagram for YBa 2 Cu 306.85 that extends out to both a and 
P — 45°, and 10 T in applied field. Strongly contrasting hqH- 
dependent behavior is observed for the two explored rotation 
directions. In general increasing a enhances the field stabil¬ 
ity range of structure B. Furthermore, for sufficiently large 
a the jUo//-dependent behavior evidences two complete first- 
order transitions that separate structure B from structure A at 
both lower and higher field. Increasing p leads to very dif¬ 
ferent behavior, since any sign of structure B is quickly sup¬ 
pressed, and completely gives way to structure A at all fields 
bvB> 17.5°. 


The phase diagram shown in Fig. 5 was constructed by 
tracking the /ro7/-dependence of the opening angles of both 
structures A and B for the range of explored a and p an¬ 
gles. These data are presented in Fig. 6. The overall ob¬ 
served behavior is generally consistent with the expectation 
of anisotropic London theory illustrated in Fig. 3. For struc¬ 
ture A [Fig. 6(a)], by increasing p, the value of at a par¬ 
ticular field increases, which is consistent with an increase of 
the effective superconducting anisotropy in the plane orthog¬ 
onal to the applied field. At low field, increasing a leads to a 
reduction in (p as the VL tends towards a more isotropic distor¬ 
tion. Fig. 6(b) shows that angle p of structure B also generally 
obeys the expected behavior as a consequence of the sample 
rotation. 

Interpreting these angle-dependent measurements in terms 
of effective mass anisotropy incorporated into local London 
theory provides at least a qualitative understanding of why the 
VL opening angle, and hence distortion, changes as a conse¬ 
quence of rotating the crystal with respect to pqH. Nonethe¬ 
less, similarly as for pqH || c, the most important aspects of 
the observations, the /io77-dependence of both cp and p, and 
the transitions observed between structures A and B, require 
a theoretical explanation that goes beyond the local approxi¬ 
mation inherent to the London theory. 

4. Discussion 

The purpose of this study is to shed new light on the 
sequence of poH-diiven and first-order VL structure transi¬ 
tions observed at low temperature in close to optimally-doped 
YBa 2 Cu 307 _^ with poH || c. Using SANS we have stud¬ 
ied the VL in a detwinned and lightly underdoped sample of 
YBa 2 Cu 306 . 85 - By underdoping we aimed to controllably al¬ 
ter the electronic structure compared with the case of lightly 
overdoped YBa 2 Cu 307 , both in terms of the precise Fermi 
surface morphology and the superconducting density of states 
(DOS). In theory, tuning these properties can lead to a varia¬ 
tion in the //oH-dependence of the VL structure between the 
two comnounds. 
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First we establish the experimental observations. For 
II c and at 7=2 K the following sequence of VL structure 
transitions has been reported in detwinned YBa 2 Cu 307 :^^’^^ 

jUo//ci < A < 2.3(2) T < B < 6.7(2) T < A 

where A/B denote VL structure types dehned in the same way 
as in Sec. 3. This /io77-dependent behavior is analogous to that 
reported here for YBa2Cu306.85 which for || c and 7=2 K 
is: 

yUo//ci < TP < 2.37(2) T < TP + B < 5.37(2) T < A. 

The principal differences between the two sequences concern 
the low and intermediate held phases. In the lowest held phase 
of YBa 2 Cu 306 . 85 , measurements with applied at an angle 
to the c-axis, allow us to deduce that structure TP arises due 
to residual twin-plane pinning in the sample. In a truly twin- 
free sample, structure A is expected to be present instead of 
structure TP, with this being the same low held VL structure- 
type seen in YBa 2 Cu 307 . 

In the intermediate held phase, both YBa2Cu306.85 and 
YBa 2 Cu 307 display structure B, yet only in the latter com¬ 
pound does structure B occupy the entire sample for || c. 
In YBa2Cu306.85, it remains unclear precisely why structures 
B and TP, which are spatially-separated in the sample, dis¬ 
play such a broad range of held-coexistence. Since measure¬ 
ments at hnite ajfi show that the co-existence of structures 
TP and B becomes replaced by a co-existence between struc¬ 
tures A and B, this allows us to rule out an obvious inhu- 
ence of the twin-plane pinning. Instead, the co-existence may 
rehect an extreme sensitivity of the VL structure to local 
variations in stoichiometry and homogeneity within the in¬ 
dividual crystals of the sample. Whatever the true explana¬ 
tion, the general observation is that in terms of absolute held, 
the propensity for structure B falls with underdoping, since it 
occupies a narrower held range in YBa 2 Cu 306.85 compared 
with YBa 2 Cu 307 . Future measurements on further under¬ 
doped samples may reveal a complete suppression of struc¬ 
ture B for ah || c. 

It is also informative to re-write the transition sequences in 
terms of magnetic helds normalised to values of Hc 2 {T - 0), 
quantitative estimates of which have recently become avail¬ 
able.For YBa 2 Cu 307 , Hc 2 is estimated to be 150(20) T, and 
the transition sequence written in terms of normalized held is: 

yUo//ci < A < 0.015(2) < B < 0.044(6) < A 

For YBa2Cu306.85, we use the value of //c2=70(10) T re¬ 
ported for YBa 2 Cu 306 . 86 , and the transition sequence is: 

/toTfci < TP < 0.034(5) < TP + B < 0.077(11) < A. 

This shows both VL structure transitions to occur at a higher 
fraction of Hc 2 in the more underdoped compound, even 
though the absolute transition helds are similar. We also point 
out that while our experiments made use of the highest mag¬ 
netic helds presently available for SANS, these helds extend 
only up to limited fractions of Hc 2 \ for YBa 2 Cu 307 the high¬ 
est applied held of 16.7 T corresponds to ~0.11 Hc 2 ?'^ while 
for YBa2Cu306.85 the highest applied held of 16 T corre¬ 
sponds to ~0.23 Hc 2 - From this viewpoint, while the obser¬ 
vations in the two compounds are comparable to one another 
on an absolute held scale, the measurements on the more un¬ 


ducting phase diagram. 

Next we discuss the physical origin of the VL structure 
transitions seen in both compounds for ^qH || c. The vast 
majority of theoretical work on VL structure transitions has 
been developed on model systems with perpendicular 
to a fourfold symmetry plane, since this provides the sim¬ 
plest system in which to treat the expected effects of either 
(or both) of a Fermi surface anisotropy and (7-wave nodal gap 

symmetry. A 1^7 p^j. perpendicular to a twofold sym¬ 
metry plane as appropriate for orthorhombic YBa 2 Cu 307 _ 3 , 
the available theory is limited^^^AS f^j- ^j-g aware, 

completely non-existent for the case where is perpendic¬ 
ular to a twofold plane that includes a nodal gap anisotropy. 
We further mention that, since held-driven and still unex¬ 
plained VL transitions have also been observed in other un¬ 
conventional superconductors for helds applied perpendicular 
to a twofold symmetry plane,^®’^° understanding the observed 
VL transitions for such a symmetry setting assumes an added 
importance from the more general perspective. In what fol¬ 
lows we interpret the observations in YBa 2 Cu 307_5 in terms 
of the expectations of available theory. 

As mentioned previously, local London theory provides a 
useful starting point for understanding the coupling between 
the VL and the crystal anisotropy in high-A: superconduc¬ 
tors.In this approach, an effective mass anisotropy is in¬ 
corporated into the free energy via a second rank tensor 
The inclusion of j can lead to distortions of isotropic hexag¬ 
onal VL structures at no energy cost, and thus local theory 
can provide an explanation for the distortions of the hexago¬ 
nal VL domains that we observe at low held in YBa 2 Cu 306 . 85 - 
However, local theory has well-noted shortcomings since no 
preferred VL orientation is expected for ^qH || c, and con¬ 
sequently no yUo7/-induced transitions in the VL coordination 
can be expected. Since in all real materials the VL displays 
a preferred orientation at ah helds, it is necessary to consider 
higher-order effects beyond local theory, i.e. nonlocal effects, 
to explain the observations. 

The inclusion of nonlocal corrections in VL theory pro¬ 
vides a means for the host system anisotropies, such as the 
Fermi velocity and superconducting gap, to couple to the VL 
properties. Consequently, nonlocal theory can be used to ex¬ 
plain both observed preferred VL orientations, and ^iqH- and 
temperature-driven VL structure transitions. Extensive theo¬ 
retical work has been conducted exploring the role of non- 
locality on the VL structure that involves varying degrees of 
approximation. Common approaches have involved develop¬ 
ing nonlocal corrections or higher-order gradients into ana¬ 
lytic phenomenological London and Ginzburg-Landau mod¬ 
els, respectively. 47 most desirable theoretical ap¬ 

proach for calculating the VL properties are hrst-principle 
calculations within the framework of quasiclassical Eilen- 
berger theory. 51 jjj approach nonlocality is catered 
for by default, and the Eermi surface and superconducting 
gap anisotropies can be treated in equal measure. Eurther- 
more, calculations can be done that are valid at any held 
or temperature. Erom such calculations it is well-established 
that for perpendicular to a four-fold symmetry plane, 
anisotropies arising from both the Eermi velocity and a nodal 
superconducting gap can determine both the VL orientation, 
and explain yUo/Z-driven VL structure transitions.^^’ 
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culations reported in Ref.^® compared the expected iiqH- 
dependence of the VL structure for cases where only a 
four-fold Fermi surface anisotropy, or only a ^f-wave gap 
anisotropy influences the VL structure. The calculations show 
that a ‘generic’ low //o^f-dependent behavior of the VL struc¬ 
ture can be expected regardless of the origin of the anisotropy. 
Briefly, the expectation is that with increasing held a first- 
order transition between a low held hexagonal VL and a 45° 
rotated hexagonal VL occurs at intermediate field. Further in¬ 
crease of the held causes the intermediate held VL to dis¬ 
tort smoothly before a second-order lock-in transition occurs 
into a square VL. The orientation of the square VL is de¬ 
termined by the symmetry of the anisotropy included in the 
model, with nearest neighbours aligned along the direction 
of the Fermi velocity minima, or the nodal direction of the 
gap. This generic expectation has good support from exper¬ 
iments, and can describe the evolution of the low held VL 
in VaSi,^^ borocarbides (Er,Y,Lu)Ni 2 B 2 C systems,and 
heavy-fermion CeCoIns.^^’^^ 

From the above discussion, it is clear that when the direc¬ 
tions of the Fermi velocity minima and gap nodes are close 
or coincide, attributing either source of anisotropy as con¬ 
trolling the VL orientation is challenging; just such a situa¬ 
tion arises for the case of two-fold symmetric YBa 2 Cu 306.85 
and YBa 2 Cu 307 . In Fig. 7 we show the Fermi surfaces of 
each compound calculated within the local-density approxi¬ 
mation (LDA) that is strictly valid only for zero-field. Each 
Eermi surface is dominated by large quasi-cylindrical bond¬ 
ing and antibonding bands arising from the Cu 02 planes, and 
a quasi one-dimensional CuO chain band dispersing along 
F-X that touches the anti-bonding band. The anisotropy of 
the dominant bonding and antibonding bands in the Bril- 
louin zone is expected to yield Fermi velocity minima along 
(110) directions. However, a contribution due to supercon¬ 
ducting chain states will shift the zone-averaged Fermi ve¬ 
locity minima towards (100). Concerning the gap symmetry, 
both phase-sensitive^® and tunnelling®® measurements done 
on YBa 2 Cu 307_3 show that the zero-field superconducting 
gap is of dj, 2 _y 2 + s symmetry. The d.^ 2 ^y 2 component with 
nodes along (110) dominates, yet the finite i-wave component 
induced by the crystal orthorhombicity causes an observed 
shift in the nodal positions also towards (100). Therefore, by 
itself each of Fermi surface and gap anisotropy are expected 
to dictate a similar field-dependence of the VL structure in 
YBa2Cu307_5. 

Bearing in mind the expected similar influences of the 
Fermi surface and gap anisotropies on the /^o7/-dependent VL 
structure, it is difficult to reconcile the observation of succes¬ 
sive first-order reorientation transitions like those observed in 
YBa 2 Cu 307_3 with the expectation of any theory. As men¬ 
tioned above, according to the calculations presented in Ref.^® 
each type of anisotropy by itself could cause the expected 
generic field-evolution of the VL structure that displays a sin¬ 
gle low held first-order reorientation transition, and a higher 
held second-order lock-in transition. In the case of a two-fold 
symmetry in the Fermi velocity, a similar generic sequence 
was found from a nonlocal London theory,although it was 
shown later that the high-field continuous transition expected 
within the model must become forbidden on grounds of sym¬ 
metry and in reality likely becomes replaced by a smooth 

'y'i, Til t • I—'* r* r-i~» y—^ y~v 



Fig. 7. (Color) In-plane Fermi surfaces of (a) YBa2Cu306.85 and (b) 
YBa 2 Cu 307 obtained from standard LDA calculations using the DMoP 
code.^^“®'^ The standard settings included DND numerical local orbital ba¬ 
sis sets, and a 4 4 2 F centered A:-points mesh. The Fermi surface contours 
were evaluated on a 50 x 50 mesh in the quarter plane. The calculations 
for the 06.85 stoichiometry was simulated in the YBa 2 Cu 307 cell, but with 
a partial band-filling corresponding to 0.3 holes per unit cell, (c) and (d) 
show the calculated Fermi surface cuts in the plane peipendicular to 
for YBa2Cu306.85T and for the cases of a=A5° andy0=45° respectively. In all 
panels, direction F-X is parallel to the crystal a* -axis, while F-Y is parallel 
to the crystal b*-axis. The black lines denote Fermi energy contours, while 
the blue (red) lines denote energy contours at Ep-\ mHa (Ep+l mHa), where 
1 mHa is 0.001 Hartree (Ha) units. 


the field-dependence of the opening angle for structure A, <p 
varies smoothly over the high-field range, and above 90° tends 
toward a saturation orientation at the very highest fields. 

On considering the above discussion, here we propose that 
the high field first-order transition seen in YBa 2 Cu 307_5 cor¬ 
responds to the “generic” first-order transition expected in 
theory that can in principle be driven by either the Fermi 
surface or gap anisotropy."® In contrast, we suggest that the 
low field transition arises instead as a consequence of a jUq//- 
dependence of the in-plane electronic structure particular to 
YBa 2 Cu 307 _ 3 , namely the superconducting DOS associated 
with the CuO chains. 

Theoretical work aimed at exploring the role of the chains 
in the vortex state of YBa 2 Cu 307_3 included calculations of 
the magnetic-field dependence of the coherence length, and 
hence vortex core size, in a bilayer model composed of plane 
and chain layers.®‘^ The calculations revealed that at very low 
field, the vortex core shape in the chain layer is orthorhombi- 
cally distorted so that it is two times longer along the chain 
direction (b-axis) than the a-axis. On the other hand, there is 
little orthorhombicity of the vortex core shape in the Cu 02 
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plane. Nonetheless, the net core shape is significantly longer 
along the b-axis than the a-axis. By increasing the field, both 
the core size in the chains, and hence the net core anisotropy, 
are expected to reduce with increasing field up to a critical 
field denoted B*. The authors estimated B* ~ 1.5 T as denot¬ 
ing a critical field above which the vortex cores in the chains 
overlap. Physically this amounts to the quenching of a compo¬ 
nent of a small superconducting gap on the chains, and hence 
a modification in the superconducting DOS. We note that B* is 
in reasonable agreement with the low VL structure transition 
field seen in both YBa 2 Cu 306.85 and YBa 2 Cu 307 , leading us 
to expect that this transition is driven by a //o//-dependent 
suppression of the core anisotropy. Verification of this expec¬ 
tation requires careful numerical calculations that make use of 
realistic material parameters for YBa 2 Cu 307 -a. In particular, 
calculating the spatial distribution of the local DOS near the 
Fermi energy, and its juof^^-dependence, are key for exploring 
the stability of the possible VL structures. On the experimen¬ 
tal side, verification could be obtained from a direct measure¬ 
ment of the low-field-dependent vortex core anisotropy by us¬ 
ing scanning tunnelling microscopy.'^’® 

The proposal that the low field transition is driven by a 
/^o^f-dependent core anisotropy can be shown as consistent 
with the SANS measurements at finite ajfi. From Fig. 5 we see 
the tendency for the low field transition to shift to lower fields 
for increasing a, and higher fields for increasing p. Qualita¬ 
tively this behavior can be understood by considering how the 
angle evolution of the Fermi surface morphology affects the 
effective core anisotropy in the plane perpendicular to the ap¬ 
plied field. To put this on a firmer footing. Figs. 7(c)-(d) show 
further LDA calculations of the Fermi surface cuts in planes 
perpendicular to the applied field for a and P-A5°. As ex¬ 
pected, at finite a and p the bonding and antibonding bands 
become elongated along the direction perpendicular to the ro¬ 
tation compared with the case a - p — 0. On the a side, 
the increase in the effective mass along the direction includ¬ 
ing the chain Fermi surface will lead to a reduction in the net 
coherence length along this direction, and hence result in a 
reduction of the core anisotropy compared with the case of 
a - p - Q. Consequently, in accordance with the expecta¬ 
tion of Ref.,® a reduced core anisotropy favors the observed 
reduction in the transition field for finite a. On the p side, 
the converse effect is observed because the effective mass in¬ 
creases in the direction perpendicular to the chains. This can 
then be expected to lead to the observed enhancement of the 
low field-stability of structure A. 

The angle-dependent measurements also shed light on the 
underlying nature of the high field transition for poH || c. 
Were the high-held transition driven solely by a nodal gap 
anisotropy, increasing a and p would be expected to have a 
largely similar effect on the transition held since each rota¬ 
tion plane includes an anti-nodal direction. In particular, by 
rotating the in-plane gap anisotropy away from the held di¬ 
rection it can be expected that a higher magnetic held, and so 
shorter intervortex distance, is required for the relevant nonlo¬ 
cal interaction to drive the structure transition. Since like the 
low transition held, the high transition held shifts in oppo¬ 
site directions for the two rotation directions, this conhrms 
that the Fermi surface morphology plays a key role in the 
high-held transition in YBa 2 Cu 307 _ 3 . It also shows that a sig- 


DOS, likely arising from persistent CuO chain superconduc¬ 
tivity, survives until the high-held region. 

Having established that the Fermi surface morphology in- 
Huences strongly the observed VL structure transitions in 
YBa 2 Cu 307_5 for poH || c, we hnally discuss the role of the 
nodal gap anisotropy. This anisotropy is generally expected 
to become increasingly important at higher helds due to the 
increasing prevalence with vortex density of the vortex lat¬ 
tice effect (VLE).^®’'^® This physical mechanism, which de¬ 
scribes the tunnelling of quasiparticles between vortices, is a 
peculiarity of the intervortex interaction in a superconductor 
with an in-plane nodal gap, and is most efficient when vortex 
nearest neighbours are aligned along the nodal directions. In 
the theory of a fourfold symmetric and purely B-wave super¬ 
conductor, the VLB is expected to lead to the stabilization of 
a perfectly square VL structure for fields above 
Note that this expected transition field as a fraction of Hc 2 
is significantly higher than the transition fields observed by 
SANS in either YBa2Cu306.85 or YBa 2 Cu 307 . This there¬ 
fore brings into question any involvement of the nodal gap 
anisotropy in driving the high field transition. 

Nonetheless, since in our SANS experiments on 
YBa2Cu306.85 we have explored the VL for applied 
fields up to ~0.23 Hc 2 , we can expect that the observed VL 
properties in the high field range are indeed sensitive to 
the gap symmetry, even if it is does not drive the high field 
transition. From this viewpoint, it is therefore interesting 
to note that, as in YBa 2 Cu 307 ,^^ at the highest fields the 
opening angle for structure A extends beyond the 90° value 
expected for a perfectly square coordination. If the very high 
field VL coordination is indeed determined by the nodal 
gap anisotropy, then the nodal directions at high field can 
be concluded to have shifted to the opposite sides of the 
(110) directions compared with the zero-field case.^®’® Full 
clarification of this intriguing high field behavior of the VL 
structure in YBa 2 Cu 307_3 would certainly benefit from both 
new theoretical work and high field experiments. 

5. Summary 

In summary, using the small-angle neutron scattering 
(SANS) technique we have studied the magnetic field- 
dependence at 2 K of the VL structure in a sample of de- 
twinned and lightly underdoped YBa2Cu306.85- For magnetic 
fields applied parallel to the crystal c-axis {poH || c), mea¬ 
surements were conducted for fields up to 16 T, and we ob¬ 
served two field-driven and first-order VL structure transi¬ 
tions at 2.37(2) T and 5.37(2) T. These transitions are anal¬ 
ogous to those reported in detwinned and lightly overdoped 
YBa 2 Cu 307 at 2.3(2) T and 6.7(2) T, and it can be expected 
that the physics behind the transitions in the two compounds 
is largely similar. In YBa2Cu306.85 further measurements at 
2 K, and up to 10 T in applied field were done after rotating 
the direction of applied field away from the c-axis. Dependent 
on whether the field is rotated around the orthogonal a*- or b*- 
axes, we observed a strong and opposing angle-dependence of 
the two transition fields, and consequently markedly different 
VL structure phase diagrams as functions of both field and 
rotation angle. 

From our observations, and in conjunction with available 
theory, we have argued here that the low-field VL structure 
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a consequence of a strong low field-dependence of a vortex 
core anisotropy that arises from the CuO chain superconduct¬ 
ing density of states. Due to the proposed involvement of the 
chain states, this low field transition can be considered as par¬ 
ticular to YBa 2 Cu 307 - 5 . 

In contrast, and once the pronounced low-field core 
anisotropy is suppressed on increasing the field, we propose 
that the high-held VL structure transition in both compounds 
more likely corresponds to the ‘generic’ hrst-order VL re¬ 
orientation transition expected in theory to be driven either 
by an anisotropy in the Fermi surface or a nodal gap. Al¬ 
though we are unable to unambiguously assign a particular 
anisotropy as driving the high held transition, the contrasting 
angle-dependences of the transition helds on the direction of 
held-rotation away from the c-axis allow us to demonstrate 
that this transition is indeed dependent on the Fermi surface 
morphology. By again comparing our experimental results to 
those of theory, we expect that the clearest evidence for an 
inhuence of the nodal gap structure on the VL properties is 
most likely away from the second transition held, and instead 
more in the range of the highest accessible helds. This can be 
checked by further held-angle SANS VL measurements done 
for very high helds above 10 T. 

Finally, and perhaps most strikingly, on the approach to the 
highest explored helds in YBa2Cu306.85 with ^qH || c= 16 T, 
the precise held-evolution of the VL coordination is ob¬ 
served to tend towards a saturation orientation. This very high 
held VL coordination is not consistent with being stabilized 
by either the established superconducting gap or Fermi sur¬ 
face anisotropies at zero-held. This ultimately demonstrates 
the strong magnetic held-dependence of the superconduct¬ 
ing density of states in YBa 2 Cu 307_5 that is rehected by the 
VL structure phase diagram. Overall it can be stated that the 
accumulated SANS observations on YBa 2 Cu 307_(5 illustrate 
the richness of the VLs in even ‘prototypical’ unconventional 
superconductors, but which still call for further experimen¬ 
tal and theoretical studies aimed at clarifying the mechanisms 
that dictate the novel held-dependent properties. 
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